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Abstract 

We propose a semi-empirical nuclear mass formula based on the macroscopic-microscopic method 
in which the isospin and mass dependence of model parameters are investigated with the Skyrme 
energy density functional. The number of model parameters is considerably reduced compared 
with the finite range droplet model. The rms deviation with respect to 2149 measured nuclear 
masses is reduced by 21%, falls to 0.516 MeV. The new magic number = 16 in light neutron-rich 
nuclei and the shape coexistence phenomena for some nuclei have been examined with the model. 
The shell corrections of super-heavy nuclei are also predicted. 
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I. INTRODUCTION 



The nuclear mass is of great importance not only for various aspects of nuclear physics, 
but also for weak-interaction studies and astrophysics [1]. In nuclear physics, it is helpful to 
study the nuclear symmetry energy and the synthesis of super-heavy nuclei by considering 
the more than 2000 measured nuclear masses. Theoretically, the mass of an atomic nucleus 
can be calculated by the macroscopic-microscopic method (such as the finite-range droplet 
model [2]) or the microscopic approaches (such as the Hartree-Fock Bogoliubov approach 
[3, 4]) or some other mass formulas [5]. The best mass formulas at present can reach about 
0.6 MeV in the rms deviation for the usual data set of 2149 measured masses of nuclei (N 
and Z > 8) [6] with about 24 ~ 30 model parameters. Compared with the microscopic 
Hartree-Fock (HF) approaches, the macro-micro model is much faster in the calculation of 
the nuclear masses for the whole nuclear chart which provides a possibility for performing a 
large scale nuclear mass calculations to refine the model parameters and to explore the global 
behavior of nuclei. However, there are two crucial points in the macro-micro method should 
be further studied. The first one is that the consistency of the model parameters between 
the macroscopic and microscopic parts in the macro-micro method should be improved. 
It is known that although the finite-range droplet model (FRDM) is widely used in the 
calculations of nuclear mass, the parameter values in the calculation of the microscopic 
shell corrections are different from the corresponding values used in the macroscopic part 
of the model [1]. This less consistency between the macroscopic and microscopic parts 
may considerably reduce the credibility of extrapolations of the macroscopic-microscopic 
approach. On the other hand, with the great development of the experimental facilities for 
the study on super-heavy nuclei and nuclei far from the /3-stability line, the influence of 
isospin effects on the nuclear mass formula attracted a great attention and should be given 
a better consideration. Based on above discussions an improved nuclear mass formula of 
self-consistently considering the isospin effects in both macroscopic and microscopic parts 
would be necessary to be established for providing a large scale nuclear mass calculations. 

To investigate the consistency of the model parameters between the macroscopic and 
microscopic parts in the macroscopic-microscopic approach and isospin dependence of the 
model parameters, the Skyrme energy density functional approach together with the ex- 
tended Thomas-Fermi (ETF) approximation [7, 8] is used. It is known that the energy 
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density functional theory is widely used in the study of the nuclear ground state which 
provides us with a useful balance between accuracy and computation cost allowing large 
systems with a simple self-consistent manner. With the Skyrme energy density functional 
approach, we systematically investigate some ground state properties of nuclei, such as the 
nuclear symmetry energy coefficient, the deformation energy and the symmetry potential, 
which are helpful to improve the macro-micro method. Based on these calculations, we 
propose a semi-empirical nuclear mass formula by taking into account the isospin and mass 
dependent model parameters. The paper is organized as follows: In Sec. II, the proposed 
mass formula is introduced. In Sec. Ill, some calculation results are presented. Finally, a 
summary is given in Sec. IV. 

II. THE MODEL 

In this section, we first introduce the macroscopic part of the mass formula. Then, the 
infiuence of nuclear deformation on the macroscopic energy of nucleus are investigated with 
the Skyrme energy density functional approach and the single particle potential used in 
the calculation of the microscopic shell correction is introduced. In addition, the symmetry 
potential and the symmetry energy coefficient of nuclear matter is also investigated. Finally, 
the parameters adopted in the model are presented. 

A. Modified Bethe-Weizsacker Mass Formula 

We start with the macroscopic- microscopic method [2, 9]. The total energy of a nucleus 
can be calculated as a sum of the liquid-drop energy and the Strutinsky shell correction AE", 

E{A, Z, f3) = Eld(A, ^) n (1 + ^^^fc') + ^' 

k>2 

The liquid drop energy of a spherical nucleus Ei,d{A, Z) is described by a modified Bethe- 
Weizsacker mass formula [10], 

Eld(^, Z) = a, A + a,A^'^ + a,^i|— H (l - Z'^/s) + a.y^I^A + a^airA-^'^Snp (2) 
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with isospin asymmetry / = (A^ — Z)/A. The pairing term proposed in [11] is adopted, with 

2 — |/| : and Z even 

|/| : N and Z odd 

1 - |/| : N even, Z odd, and iV > Z 

1 - |/| : N odd, Z even, and < Z 

1 : even, Z odd, and N < Z 

1 : A^ odd, Z even, and N > Z 



(3) 



In this work, the symmetry energy coefficient of finite nuclei is written as, 

1 
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(4) 



based on the conventional surface-symmetry term [12, 13] of liquid drop model, with a small 
correction term for description of isospin dependence of a^ym- The sensitive dependence of 
symmetry energy coefficient on the asymmetry of nucleus, especially that agym increases with 
increasing proton fraction of the system is also found in [14]. The introduced I correction 
term approximately describes the Wigner effect [2] of heavy nuclei. For a heavy nucleus 
near the /3-stability line {A\I\ ^ 2 ^ |/|), the introduced / term in asym roughly leads to a 
correction oc |/| (known as Wigner term) to the binding energy of the nucleus. Compared 
with the case without the / term being taken into account, the rms deviation of nuclear 
masses defined as = ^ ^Mixp — ) from the measured masses AME2003 [6] for 
the 2149 nuclei (A^ and Z > 8) is reduced by about 6%. Furthermore, we find when 
the isospin dependence of symmetry energy coefficient is taken into account, the obtained 
optimal Cgym increases about 3 MeV and up to 29 MeV which is close to the symmetry 
energy coefficient of nuclear matter at saturation density obtained from the Skyrme energy 
density functional. 

The Coulomb exchange correction and surface diffuseness correction to the Coulomb 
energy is approximately taken into account as the term Z^"^!^ . In addition, the terms 6^ in 
Eq.(l) which are obtained according to the Skyrme energy density functional (the detailed 
discussion is in the next subsection) describe the contribution of nuclear deformation to the 
macroscopic energy. 
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FIG. 1: (Color online) (a) Energy of ^^O, ^®Ca and 208p]-, ^j^]-^ respect to /32 deformation. Here, the 
values of Eq are negative. The circles and the solid curves denote the results of SkM* interaction 
and of a formula E/Eq = 1 + ^2/32) respectively, (b) The value of 62 obtained with SkM* as a 
function of mass number. 

B. Influence of Nuclear Deformation on the Macroscopic Energy 

For the deformation of nuclei, we only consider axially-deformed cases. In this work, only 
(^2 and Pi deformations of nuclei are taken into account. We first investigated the energy of 
a nucleus with respect to a (32 deformation based on the Skyrme energy density functional 
together with the extended Thomas- Fermi approximation (ETF) [7, 8]. The procedure is 
as follows: The total energy of a nucleus can be expressed as the integral over the Skyrme 
energy density functional 'H(r) [15]. Given a density functional p(r), one can calculate the 
corresponding energy via E = J T-L[p[r)]dr under the ETF approximation. We first obtain 
the binding energy Eq and the spherical Woods-Saxon density distributions of a nucleus 
with the approach in [8]. Then, with the same procedure we calculate the energy E{P2) of 
the nucleus with a quadrupole deformed Woods-Saxon density distribution of the nucleus in 
which the central density and the surface diffuseness remained unchanged. Fig. 1(a) shows 
the calculated energy of ^^O, ^®Ca and ^osp]^ ^ function of P2 deformation with the 
SkM* interaction [15] (denoted by circles). The solid curves denote the results of a formula 
E/Eq = 1 + 62/31 which the value of 62 is obtained by fitting the open circles. One can 
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FIG. 2: (Color online) (a) Energy of ^'^^Pb with respect to /32 and /34 deformation with Skyrme 
energy density functional approach, (b) The value of bk as a function of mass number. The 
scattered symbols denote the obtained curvatures of the parabolas with the Skyrme force SkM* for 
a number of nuclei. The solid curves denote the corresponding results with an empirical formula 
(6). 

see that the parabola approximation to the change of energy with P2 is acceptable. Fig. 1(b) 
shows the value of 62 as a function of the mass number. The crosses denote the results of 
the SkM* interaction for a number of nuclei along the /3-stability line. We find that the 
dependence of 62 on the mass number A can be reasonably well described by a formula 

b2 = giA'/' + g2A~'/\ (5) 

This form of mass dependence of 62 is therefore adopted in the proposed mass formula and 
the optimal values of gi and g2 are finally determined by the 2149 measured nuclear masses 
[6]. 
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To take into account the influence of the higher-multipole deformation of nuclei, we inves- 
tigate the change of energy of a nucleus with respect to a certain set of nuclear deformation 
parameters with the Skyrme energy density functional approach mentioned previously. In 
Fig. 2 (a), we show the energy of ^osp]^ function of P2 and with the Skyrme force 
SkM*. We flnd that the influence of nuclear P4 deformation on the nuclear energy can be 
roughly described by a parabola at small deformations. For other higher-multipole defor- 
mations, the parabola approximation can also be applied for small deformation cases, and 
the proposed model can be easily extended to consider other higher-multipole deformations. 
Furthermore, we notice that the curvature of the parabola for a given (3^ deformation can 
be approximately described by an empirical formula 

6.= (|)^iAV^+Q)^.A-V3, (6) 

which is an extension of the formula (5). In Fig. 2(b), we show the mass dependent curvatures 
of the parabolas. The crosses, the open circles, the solid squares and the triangles denote 
the obtained curvatures 62? b^, 64 and feg of the parabolas with the Skyrme energy density 
functional approach for a number of nuclei along the /3-stability line, respectively. The solid 
curves denote the corresponding results with the empirical formula (6) taking gi = 0.0074 
and g2 = —0.38. One can see that the curvatures of the parabolas can be reasonably well 
described by the empirical formula which greatly reduces the computation time for the 
calculation of deformed nuclei. 

C. Single-particle Potential in the Microscopic Part 

In the microscopic part, the shell correction 

AE = ciE,h, (7) 

is obtained by the traditional Strutinsky procedure [16] by setting the smoothing parameter 
7 = 1.2hwQ and the order p = 6 of the Gauss-Hermite polynomials. Where, Esh = -Esh(-P) + 
Esh{N) i.e. sum of the shell energies of protons and neutrons, hwo = 41yl~^/'^MeV is the 
mean distance between the gross-shells. In this work, we introduce a scale factor ci to the 
shell correction. This additional parameter is used to adjust the division of the binding 
energy between the macroscopic part and the remaining microscopic correction. It is known 
that a similar scale factor is usually introduced to the liquid-drop part [17] or the shell 
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correction part [18] to adjust the division between the two parts for giving better results in 
the calculation of fission barrier. It is necessary to investigate the influence of this parameter 
on the nuclear masses. We find that the rms deviation for the 2149 nuclear masses can be 
somewhat reduced with the introduced factor ci compared with the case setting ci = 1. 

To obtain the shell correction Ai?, we execute a computer code WSBETA [19] to calculate 
the single particle levels of an axially deformed Woods-Saxon potential and then perform 
the Strutinsky procedure. The single-particle Hamiltonian in the code WSBETA is written 
as 

H = T + V + V,,,,, (8) 

with the spin-orbit potential 



K.o. = -A(^^J VV-{axp), (9) 

where A denotes the strength of the spin-orbit potential. In this work, we set A = Aq (l + ^) 
with Ni = Z for protons and Ni = N for neutrons. Here, the isospin-dependent spin-orbit 
interaction strength is obtained based on the Skyrme energy-density functional in which the 
spin-orbit potential is usually expressed as 

= VoV(p + p,) -iaxp)^ Vo (l + ^]\/p-iaxp), (10) 



with the nucleon density p = Pp + Pn and the spin-orbit strength Wq- M in Eq.(lO) is the free 
nucleonic mass, a and p are the Pauli spin matrix and the nucleon momentum, respectively 
[19]. The central potential V is described by an axially deformed Woods-Saxon form 

1 + exp(^^) 

Where, the depth Vq of the central potential {q = p for protons and q = n for neutrons) is 
written as 

Vq = Vo± Vsl (12) 

with the plus sign for neutrons and the minus sign for protons. Vs is the isospin-asymmetric 
part of the potential depth. We assume Vg = agym in this work (detailed study of the relation 
between Vg and a^y^ is given in the following part of this section). TZ defines the distance 
from the origin of the coordinate system to the point on the nuclear surface 

n{e) = coR [1 + /32Y2o{e) + /3,Y,o{e) + ...], (13) 



with the scale factor cq which represents the effect of incompressibihty of nuclear matter in 
the nucleus and is determined by the so-called constant volume condition [19]. Yim{0, </>) are 
the spherical harmonics. R = tqA^^^ and a denote the radius and surface diffuseness of the 
single particle potential, respectively. Here, we assume and set the radius and diffuseness of 
the single particle potential of protons equal to those of neutrons for simplicity. For protons 
the Coulomb potential is additionally involved (see [19] for details). 

D. Symmetry Potential and Symmetry Energy CoefRcient 

The relation between the isospin-asymmetric part Vg of the single particle potential depth 
in the microscopic part and the symmetry energy coefficient in the macroscopic part is 
investigated based on the Skyrme energy density functional together with the ETF approach. 
In this approach, the central one-body potential is described by Vg = with the energy 
density functional e{r) (see Eq.(9) in Ref. [7] for details). The difference between the neutron 
{q = n) and proton (g = p) potentials of nuclear matter is written as 

= 2B2p6 + 2B8p''^'6 + BiCkp'^^6 + 0{6^) (14) 

with the kinetic energy density Tg which can be expressed as Tg = |(37r^)^/^pg'^^ in the 
Thomas- Fermi approximation, the isospin asymmetry 6 = (p„ — Pp) / p, and the coefficient 
Cfc = (37r^/2)^^^. B2, Bg and -B4 (notations in [7]) are some combinations of Skyrme param- 
eters, given by B2 = -|to(| + a;o), ^8 = - j^^ad+a^s) and B4, = -|[ti(| +Xi) -tad + a;2)]- 
The symmetry potential Vsym may be written as 

Vsym = ^^^^ = B2P + Bsp-^' + Ib,c,p'/' + 0{6'). (15) 

The symmetry energy coefficient of nuclear matter J is written as [20] 

,.is.p4B,p"«-le.c./'».i(^)c.,v. (16) 

with Qg = 3tiXi— ^2(4+5x2). The 6s term and the last term of Eq.(16) give the contributions 
of the effective-mass [7] and the kinetic energy to the J, respectively. From the above 
equations for Vsym and J, one can get the relation between them, 

J = \Vsym - ^[6. + 6i?4]c.p^/^ + \ (J^^ c,p"\ (17) 




J (MeV) (N.z)/A 



FIG. 3: (Color online) (a) Symmetry potential Vsym of nuclear matter with 78 Skyrme forces, 
(b) The difference between the potential depth of neutrons and of protons as a function of mass 
asymmetry. The circled curve and the triangled curve denote the results fitted to the calculated 
results with SLy4 (short red dashes) and SkM* force (small gray squares) for a large number of 
nuclei, respectively. 

A similar equation is previously proposed in [21] based on perturbation theory, 



dVojkr, 
dkr. 



1 / 



3 V2M 



ki. (18) 



Due to the uncertainty of choosing the interaction parameters, there exists a large un- 
certainty for the value of Kj/m in different models. In Fig. 3 (a), we show the calculated 
symmetry potential Vsym of nuclear matter with 78 Skyrme forces. The Vsym has a value 
of about 10 ~ 50 MeV according to the calculations. Brueckner-Hartree-Fock calculations 
show that the value of Vsym is about 25 MeV [22]. These calculations indicate that the value 
of Vsym is comparable to that of the symmetry energy coefficient J which is about 30 MeV. 

For finite nucleus, the isospin-asymmetric part Vs of the single particle potential should 
be slightly different from the value of Vsym- With the density distributions of nuclei obtained 
in [8] , we calculate the potential depth of protons and neutrons for a large number of nuclei. 
We find that the difference Vn — Vp increases linearly with the isospin asymmetry I (see 
Fig.3(b)). The average value for the isospin-asymmetric part Vs can be obtained by linearly 
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TABLE I: Model parameters of the mass formula. 



parameter 


WS 


a^, (MeV) 


-15.5841 


as (MeV) 


18.2359 


ac (MeV) 


0.7173 


CsymiMeY) 


29.2876 


K 


1.4492 


apair{MeV) 


-5.5108 


91 


0.00862 


92 


-0.4730 


Cl 


0.7274 


Vo (MeV) 


-47.4784 


ro (fm) 


1.3840 


a (fm) 


0.7842 


Ao 


26.3163 



fitting the calculated results. The obtained values of Vg are 17.0 and 26.9 MeV with SLy4 
[23] and SkM* [15] force, respectively. In [24], the authors found that the experimental 
Fermi energies of a number of magic nuclei can be well described with a value of 23.2 MeV 
for Vs. The asym in this work has a value of about 23 ~ 24 MeV for heavy nuclei which 
is roughly comparable to the obtained values of In the first round of searching for the 
optimal parameters of the proposed mass formula, we treat K as a free parameter and find 
that the obtained value of Vs is very close to that of asym- So we empirically set and assume 
Vs ~ CLsym in the improved mass formula for simplification. 

E. Model Parameters 

From the above discussions, one can see that the macroscopic and microscopic parts in 
the proposed mass formula are closely connected to each other through the coefficient asym of 
the symmetry energy and other isospin dependent model parameters. The number of model 
parameters is considerably reduced compared with the finite range droplet model (FRDM) in 
which the number of parameters is about 31 [1]. Here, we have 13 independent parameters a„. 
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TABLE II: rms a deviations between data AME2003 [6] and predictions of several models (in 
MeV). The line a{M) refers to all the 2149 measured masses, the line cr(5„) to the 1988 measured 
neutron separation energies 5^. The calculated masses with FRDM are taken from [2]. The masses 
with HFB-14 and HFB-17 are taken from [3] and [4], respectively. 





FRDM 


HFB-14 HFB-17 


WS 


a{M) 


0.656 


0.729 


0.581 


0.516 


o"(5'„) 


0.399 


0.598 


0.506 


0.346 


Np 


31 


24 


24 


15 



Cpair, 5'!, 5'2, Ci, Vo, ro, a, Ao for the nuclear mass. By varying these parameters 
and searching for the minimal deviation of the 2149 nuclear masses from the experimental 
data, we obtain a parameter set labeled as WS which is listed in Table 1. To find the 
minimal energy E{A,Z,/3) with respect to a set of deformation parameters for a given 
nucleus, the downhill searching method is adopted. We re-execute the downhill algorithm 
for several times starting from different initial deformation parameters in order to find the 
lowest energy of a nucleus from some possible local minima on the energy surface E{A, Z, 
In the calculation of nuclear masses with the obtained binding energies, the electron binding 
energies are not included. In the parameter searching procedure, the downhill searching 
method and the simulated annealing algorithm [25] are incorporated. The former is used 
for the parameters of the microscopic part, while the latter is for the macroscopic part. 

III. RESULTS AND DISCUSSION 

In this section, we first show the calculated rms deviations of the nuclear masses and of 
the neutron separation energies. In addition, the change of magic number in light neutron- 
rich nuclei and the shape coexistence phenomena for some nuclei have been checked with 
the model. Then, the shell corrections of super-heavy nuclei and the location of the center 
area of the super-heavy island are investigated with the proposed mass formula. 

A. Test of the Model 

The corresponding rms deviations of nuclear masses for the 2149 measured nuclei with 
the parameter set WS is listed in Table 2. In addition, the results of FRDM and Hartree- 
Fock Bogoliubov (HFB-14 [3] and HFB-17 [4]) are also listed for comparison. Np denotes 

12 




FIG. 4: (Color online) (a) Deviations between the calculated nuclear masses from the experimental 
data, (b) Calculated shell corrections AE of nuclei (crosses). The squares denote the microscopic 
energy of nuclei with the FRDM model (column Emic of the table of Ref . [2] ) . 

the corresponding number of parameters used in each model. Compared with the FRDM, 
the rms error for the 2149 nuclear masses is considerably reduced with WS, from 0.656 to 
0.516 MeV. The number of parameters in the model is reduced from 31 to 15 (including the 
two parameters 7 and p used in the Strutinsky procedure). One should note that several 
(about 12) of the FRDM parameters were prefixed by considerations other than mass-like 
data before making the fit [1, 2] and also that the fit in the FRDM included data on fission 
barriers in addition to masses. In this work, only the precisely measured nuclear masses are 
involved in the fit. Compared with the standard Hartree-Fock Bogoliubov (HFB) approach, 
the CPU time used in the calculation of nuclear mass table is much shorter with the proposed 
mass formula. The obtained rms error for the 1988 measured neutron separation energies 
Sn with our model is obviously smaller than those of HFB calculations [3, 4]. 

Fig. 4 (a) shows the deviations between the calculated nuclear masses in this work from 
the experimental data. In Fig. 4(b), we show the calculated shell corrections AE of nuclei 
with our model and the microscopic energy (mainly including the shell correction and the 
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FIG. 5: (Color online) Potential energy surface E{(3\ 



i2,/34) of 82Sr. 



deformation energy) obtained in the finite-range droplet model. For intermediate and known 
heavy nuclei, the results of the two approaches are comparable and both of them reproduce 
the known magic numbers very well. The deviations are large for light nuclei and super- 
heavy nuclei. Our calculations show that the shell corrections of nuclei with about = 16 
are much larger (in absolute value) than those from the FRDM. As an example, the shell 
correction of ^^O is calculated and has a value of —4.6 MeV with WS. Experimentally, it 
is thought that ^^O is a doubly magic nucleus from the observed decay energy spectrum 
and the high-lying first excited 2+ state (above 4.7 MeV) [26], which is consistent with our 
calculations. The obtained shell corrections with WS for ^°C, ^^C [27] and ^^N are —4.2, 
—5.1 and —4.3 MeV, respectively. Some theoretical and empirical studies [26, 28] have shown 
that in the neutron-rich nuclei the magic numbers such as = 14 or 16 can arise, which is 
in agreement with our calculations. It is known that the shell correction strongly depends 
on the single particle potential adopted. The isotopic dependence of the spin-orbit strength 
and the symmetry potential adopted in this work is different from that in the FRDM, which 
leads to the different shell correction from the two models. Our results for the neutron-rich 
nuclei with about A^ = 16 look more reasonable qualitatively. 

To further test the model, we study the potential energy surface £'(/32, P^) of some nuclei. 
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TABLE III: Shell corrections of some nuclei (in MeV). The data of FRDM are taken from the 
microscopic energies Emic of the table of Ref. [2] . 





160 


240 


40Ca 


48Ca 


902r 


i32Sn 


208p]^ 


270Hs 


288]^]^4 


292114 


298114 


294116 


ws 


-0.7 


-4.6 


2.0 


-1.2 


-1.3 


-9.8 


-11.0 


-6.4 


-5.3 


-6.1 


-6.1 


-6.2 


FRDM 


2.1 


0.3 


2.3 


0.1 


-1.6 


-11.6 


-12.8 


-6.5 


-7.8 


-8.9 


-7.6 


-8.7 



In [29], the authors observed the shape coexistence phenomena for nuclei ^"^St and Kr iso- 
topes from the low and high spin states. The shape coexistence phenomena of these nuclei 
could be observed from the corresponding potential energy surface. In Fig. 5, we show the 
calculated potential energy surface of ^"^Si. The coexistence of oblate and prolate deformed 
configurations can be clearly observed. The similar coexistence phenomena for Kr isotopes 
can also be observed with our model. 

B. Shell Corrections of Super- heavy Nuclei 

The precise calculation for the shell corrections of super-heavy nuclei is of great impor- 
tance for the synthesis of new super-heavy nuclei, especially for the prediction of the location 
of the super-heavy island. Furthermore, the fission barriers of super-heavy nuclei are roughly 
estimated by the values of the corresponding shell corrections [18, 30] since the macroscopic 
fission barriers disappear at super-heavy region in general. It is known that the fission bar- 
rier is a very sensitive parameter in the realistic calculations for the survival probabilities of 
the produced compound nuclei. It is therefore necessary to investigate the shell corrections 
of super-heavy nuclei. 

In Table 3, we list the calculated shell corrections AE of some nuclei. The corresponding 
microscopic energy obtained in the FRDM are also listed for comparison. For super-heavy 
nuclei such as nucleus 292114^ ^]-^g rnicroscopic energies obtained with the FRDM are much 
lower (absolute value larger about 2 ~ 3 MeV) than our calculated AE. Because these nuclei 
are (nearly) spherical in shape according to the calculations. It follows that the deviations of 
shell energies between the two models are about 2 ~ 3 MeV for nuclei around 292114. Because 
the shell correction can not be measured directly, it is still difficult to quantitatively compare 
the reliability of model through the calculated shell corrections of nuclei. In addition, we 
study the central area of the super-heavy island based on the calculated shell energies. Fig. 6 
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FIG. 6: (Color online) (a) Shell correction energies AE of nuclei. The black squares denote the 
nuclei with microscopic energies of — (6 ~ 7) MeV in the FRDM calculations. The straight line 
passes through the areas with the known heavy magic nuclei, (b) Shell correction energies of nuclei 
in super-heavy region. The crosses denote the nearly spherical nuclei (calculated |/32| < 0.01) and 
the triangles denote the synthesized super-heavy nuclei in the "hot" fusion reactions [32, 33]. 

shows the contour plot of the calculated shell correction energies of heavy nuclei. The black 
squares in sub-figure (a) denote the nuclei with microscopic energies of —(6 ~ 7) MeV in 
the FRDM calculations. One can see that both of models give similar magic numbers for 
heavy nuclei. Fig.6.(b) shows the shell corrections of nuclei in the super-heavy region. The 
crosses denote the calculated nearly spherical nuclei (|/32| < 0.01) with WS. The predicted 
super- heavy island according to the obtained shell corrections of nuclei looks fiat. Along 
nuclei with Z = 114 or = 178, one can see a slightly deeper valley in the contour plot 
of shell corrections. The calculated deformations of nuclei demonstrate that the nuclei with 
= 184 are (nearly) spherical in shape. However, the maximum shell correction occurs at 
A^ = 178 instead of A^ = 184, which is consistent with the results in [30, 31]. The analysis 
about the shift of the shell correction from A^ = 184 to A^ = 178 is given in Ref. [30]. 
According to the calculations, the super- heavy nuclei 288,289]^]^^ produced in the "hot" fusion 
reaction ^^Ca+^^^Pu [32] (the corresponding compound nucleus is ^^^114) are close to this 
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FIG. 7: (Color online) Deviations of the calculated nuclear masses in this work from the results 
of FRDM (a) and HFB-17 (b), respectively. The calculated masses with FRDM and HFB-17 are 
taken from [2] and [4], respectively. The shades denote the region with deviations smaller than 2 
MeV. 

central area of the island. The half-lives of these nuclei are in the order of seconds [32], 
which is much shorter than those of known stable nuclei. The measured short half-lives of 
nuclei in super-heavy region seem to indicate that the shell corrections of these nuclei are 
probably not very large. 

Fig. 7 shows the deviations of the calculated nuclear masses with the proposed model from 
the results of FRDM and HFB-17. The shades denote the region with deviations smaller 
than 2 MeV. The results for highly neutron-rich heavy nuclei from the three models have 
large deviations. The results of our model are relatively close to those of HFB-17 for most 
nuclei. 

IV. SUMMARY 

In this paper we proposed a semi-empirical nuclear mass formula based on the 
macroscopic-microscopic approach. The isospin effects in both macroscopic and microscopic 
part of the formula are self-consistently considered, with which the number of model pa- 
rameters is considerably reduced compared with the finite range droplet model and the rms 
deviation of the calculated masses from the 2149 measured nuclear masses is reduced by 
21% and falls to 0.516 MeV. The CPU time used in the calculation of the nuclear masses 
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for the whole nuclear chart is much shorter than that with the microscopic mass formula 
models. At the same time the consistency of the model parameters between the macroscopic 
and microscopic parts greatly promotes the credibility of extrapolations in the macroscopic- 
microscopic approach. 

In order to extend the mass formula to super-heavy nuclei and the nuclei far from the 
/3-stability line, we pay a special attention to study the isospin and mass dependence of the 
model parameters including symmetry energy coefficient and the symmetry potential, etc. 
Those studies are based on the Skyrme energy density functional approach together with 
the extended Thomas-Fermi approximation. Since more sufficiently considering the isospin 
effects of the model parameters the formula could systematically study super-heavy nuclei 
and the nuclei far from the /3-stability line. 

To further test the model, the appearance of new magic number = 16 in light neutron- 
rich nuclei and the shape coexistence phenomena for some nuclei have been examined with 
the model. Our results are in good agreement with some experimental and theoretical 
studies. The predicted super-heavy island according to the obtained shell corrections of 
nuclei looks flat. Along nuclei with Z = 114 or A^ = 178, we flnd a relatively deeper valley 
in the contour plot of shell corrections. The shell corrections of nuclei around ^^^114 are 
about —6 MeV and much smaller (in absolute value) than the corresponding results from 
the flnite-range droplet model. The calculated nuclear masses for highly neutron-rich heavy 
nuclei from the three different models have large deviations. The results of our model are 
relatively close to those of HFB-17 for most nuclei. 
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